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The exact diagonalization and the variational cluster approximation (VCA) are used to 
study the nature of a novel Mott insulator induced by a strong spin-orbit coupling for a two- 
dimensional three-band Hubbard model consisting of the t2 g manifold of 5d orbitals. To char- 
acterize the ground state, we introduce a local Kramer's doublet which can represent a state 
with effective angular momentum J e g = \S — L\ = 1/2 as well as spin S = 1/2. Our system- 
atic study of the pseudo-spin structure factor defined by the Kramer's doublet shows that the 
JcH = 1/2 Mott insulator is smoothly connected to the 5 = 1/2 Mott insulator. Using the 
Kramer's doublet as a variational state for the VCA, we examine the one-particle excitations 
for the Mott insulating phase. These results are compared with recent experiments on S^IrCU. 

KEYWORDS: 5d electrons, spin-orbit coupling, Mott insulator, three-band Hubbard model, variational 
cluster approximation 



1. Introduction 

bd transition metal oxides have attracted much atten- 
tion because of their unique properties caused by a strong 
spin-orbit coupling (SOC) A for 5d transition elements. 
One of such materials is Sr2lr04 in a layered perovskite 
structure of K2N1F4 type. 1 ' 2 ' While Sr 2 Ir04 was first 
synthesized more than fifty years ago, 1 ' it is only in 90 's 
that its electronic properties has been studied systemat- 
ically as an analogous system to high-T c cuprate super- 
conductors. An inelastic neutron scattering experiment 
has found that the ground state of Sr2lr04 is antiferro- 
magnetically ordered with weak ferromagnetic moment, 
similar to that for the parent compounds of cuprates. 2 ' 
The early studies then concluded that the ground state 
of Sr2lrC>4 was a spin S = 1/2 Mott insulator with un- 
paired electrons occupying a half-filled d xy band, which 
is split off upward compared to d yz and d zx orbitals due 
to the crystalline electrostatic field with elongation of the 
Ir-0 bond along z direction. 3 ' However, very recently, x- 
ray scattering experiments have revealed that the ground 
state is instead close to a J e g = 1/2 Mott insulator. 4 ' 
Here, J e g = \S — L\ is an effective total angular momen- 
tum defined in the tig manifold with the orbital angular 
momentum L. Note also that the J e g = 1/2 state corre- 
sponds to the ground state in the atomic limit with large 
A. 

Motivated by these experiments, we shall study the- 
oretically the electronic and magnetic properties of 
Sr2lr04 using a three-band Hubbard model with the 
SOC. 5 ' First, we introduce a local Kramer's doublet 
which can represent a S = 1/2 state as well as a 
Jefi = 1/2 state. 6 ' Employing the exact diagonalization 
method, we show that the local Kramer's doublet can 
describe smoothly both the S = 1/2 and J c s = 1/2 
Mott insulators with varying A. This strongly indicates 
that there exists no apparent symmetry change between 

*E-mail: t-shirakawa@riken.jp 



these two extreme states. We then employ the variational 
cluster approximation 7 ' (VCA) method based on the self- 
energy functional theory 8 ' (SFT) to examine the one- 
particle excitations in the Mott insulating phase. We find 
that, for a realistic set of model parameters for Sr 2 Ir04, 
most of the unoccupied state, i.e., the upper Hubbard 
band, can be well describe by the J c s = 1/2 state. This 
result is in good qualitative agreement with the recent 
experimental observation. 4 ' We also discuss the effects of 
SOC and local Coulomb interactions on the one-particle 
excitations. 

The paper is organized as follows. After describing the 
three-band Hubbard model in Sec. 2, the local Kramer's 
doublet is introduced in Sec. 3 and compared with the 
numerically exact ground state of the model. Using VCA 
based on SFT, the one-particle excitations are studied in 
Sec. 4. The paper is concluded in Sec. 5. 

2. Model 

Since the crystalline electrostatic field is much larger 
than the SOC and the local Coulomb interactions, 10 ' the 
local electronic configuration of Ir 4+ ion in Sr2lr04 is 
the low-spin state of fog) 5 - Therefore, we consider the 
following effective three-band Hubbard model, consisting 
of three d orbitals {d xy , d yz , and d zx ), on the square 
lattice 
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Table I. Two sets of parameters used. 



Parameters 


Simplified model 


Sr2lr04 


h 


0.36 eV 


0.36 eV 


t2 


0.00 eV 


0.18 eV 


t 3 


0.00 eV 


0.09 eV 


u 


0.36 eV 


0.37 eV 


ts 


0.00 eV 


0.06 eV 


A 


0.36 eV 


-0.36 eV 


X 




0.37 eV 



+ J 'EE 

CrjSiCr^- (4) 

r a^f) 

Here, c rQ(7 (cJ, Qcr ) is the annihilation (creation) operator 
of an electron with spin a (a =^, I) and orbital a (a — 
xy, yz, and zx) at site r. Ck aa is the Fourier transform 
of c raa . is the dispersion of orbital a 

ej£ = — 2ti(cos k x + cos k y ) — t<i cos k x cos k y 

-2t 3 (cos2k x + cos 2k y ) + A, (5) 

(6) 

(7) 



-2ti cos k y — 2t§ cos k x , 



-2td cos k-r — 2t<=. cos k„ 



where t\, t%, and t 3 correspond to the hopping integrals 
for d xy orbitals located at the nearest, next nearest, and 
third nearest neighbor sites, respectively. t± and £5 are 
the nearest neighbor hopping integrals for d yz (d zx ) or- 
bital in y (x) and x (y) directions, respectively. A is an 
energy level difference between d xy orbital and the other 
orbitals (d yz and d zx ), which is naturally expected due to 
the large crystalline electrostatic field inducing tetrago- 
nal splitting. H so is the SOC term represented in the 
following matrix form 
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where s = +1 (—1) for a =f (\), and a indicates the 
opposite spin of a. Finally, we introduce, for the local 
Coulomb interactions, the intra-orbital (U) and the inter- 
orbital (U') Coulomb interactions, the Hund's coupling 
J, and the pair- hopping J', with U = U' + 2 J and J = 
J'. 9 ) The number of electrons is set to be 5 per site. 

As shown in Table I, we use two sets of parameters 
indicated by " Simplified model" and " S^IrCV . The pa- 
rameter set "Sr2ir04" is a realistic set of parameters de- 
termined by fitting the band dispersion obtained from 
the first principles calculations based on density func- 
tional theory. 10 ) 

3. Exact diagonalization study 

First, we have obtained the phase diagram for the 
simplified model with U/ti = 8 and J/U = 0.15 us- 
ing exact diagonalization technique for a 4-site cluster. 
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Fig. 1. (Color online) Structure factors C^(q) (for definition, see 
in the text) for (a) q = (0,0), (b) q = (tt,0), and (c) q = (it, it). 
(d) The optimal value of 8 for the Kramer's doublet [Eq. (9)]. 
The model parameters used are U/ti = 8, J/U = 0.15, A/fi = 1, 
and the simplified hopping listed in Table I. 



Let us briefly summarize the phase diagram of A vs. 
A. 11 ) In region (I) A/ii < -0.3 and A/ti < 0.15, d yz 
and d zx orbitals are magnetically active, whereas most 
of d xy orbital is fully occupied. The degeneracy of d yz 
and d zx orbital leads to a ferromagnetic state with anti- 
ferro orbital ordering, which is expected by the second- 
order perturbation theory from the strong coupling limit 
in the presence of Hund's coupling. 12 ' In region (II) 
-0.3 < A/ti < 0.4 and \Jt x < 0.2, the hole density 
is n xy : n yz + n zx = 1/2 : 1/2 for A = 0, and stripe-like 
spin correlations, characterized by q = (7r,0) and (0,7r), 
become dominant. Finally, in region (III) A/ii ~ 1 or 
X/ti > 0.2, an antiferromagnetic insulating phase ap- 
pears, including both S = 1/2 and J er j = 1/2 Mott insu- 
lators in the limit of A = and co, respectively. 

Let us now discuss closely the ground state proper- 
ties in region (III) . We first introduce the following local 
Kramer's doublet 

Si ""' ' ■ ' ' (9) 



a r0a = SCOS0C ixya 



V2 



'ryza 



which includes the two limiting states, namely, 5 = 1/2 
(d xy ) for 9 = and J c s — 1/2 for 9 — arctan\/2 ~ 
54.74°. In the atomic limit, the eigenstate of the highest 
level of H so in the presence of A is generally expressed 
by this state. 6 -* Here, we shall optimize 9 by minimizing 
the hole density ng a = (ipo \ a r e a o\,g a |^o)- This quantity 
should be ng a — 0.5 if a\, ga is a well-defined particle. The 
numerically obtained optimized 9 is summarized in Fig. 
1 (d) for the simplified model, where we can see that the 
optimized 9 varies smoothly with A. The deviation of ne a 
from 0.5 is found less than 1.2% for a range of A studied, 
indicating that the ground state is well described by the 
Kramer's doublet. 

Next, let us study the magnetic properties in region 
(III) of the phase diagram. For this purpose, we first 
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define the structure factor for a local angular momentum 
operator 0£f 



,iq-(i — r') 



(Vol 0» ■ 0£ |Vo) , (10) 



r,r' 



where N is the number of lattice sites, and \ipo) indicates 
the ground state. The spin structure factor Cs(q) is then 
obtained simply by setting 

' ^ (11) 



O b r = S r = 



2^ 

aaa' 



where er is the vector representation of Pauli matrices. 
The structure factor for the effective total angular mo- 
mentum Cj B{1 (q) is defined by introducing the following 
local angular momentum operator 



o> 



£ 



r T L , 



(12) 



where L is the orbital angular momentum operator. 
Note that as long as the local bases are confined within 
tig manifold, the matrix elements among ti g states are 
equivalent to those among p states apart from minus 
sign. 13 ) Finally, we define the structure factor for the 
Kramer's doublet Ckd(<z) by introducing the following 
pseudo-spin 1/2 operator 



(13) 



where a' rga is given in Eq. (9). 

The numerical results for these structure factors are 
shown in Fig. 1. For A = 0, the ground state is the 
5=1/2 Mott insulator, where the hole density of d xy 
orbital is exactly 1 whereby the model is equivalent to 
the single-band Hubbard model. With increasing A, wc 
can see in Fig. 1 (a)-(c) that the values of Gj cff (q) ap- 
proach to those of Cs{q) for A = 0. This indicates that 
the ground state in the limit of large A is represented 
simply by an antiferromagnetic ordering of J c ff = 1/2 
angular momentum. Finally, we also find in Fig. 1 (a)— 
(c) almost no A dependence on Gkd(<7)> which strongly 
suggests that the ground state for different values of A 
can be well described by a state with alternative align- 
ment of the Kramer's doublet pseudo-spin. Therefore, 
we conclude that the 5 = 1/2 and J e g = 1/2 Mott in- 
sulators, the two extreme states for A = and oo, are 
smoothly connected with no apparent symmetry change. 

4. Variational Cluster Approximation Study 

We now adopt the VCA method 7 ' based on the SFT 8 ' 
to study the low-energy one-particle excitations. This 
method takes into account precisely the effects of short- 
range static and dynamical correlations, and thus it is 
superior to a simple mean field approximation. The SFT 
introduces a reference Hamiltonian H' with the same 
two-body interactions as H but with a different one- 
body part t' ', and H' may be solved numerically ex- 
actly on a finite cluster. An approximate grand poten- 
tial for H is given in a functional form by f2(t') = 
il' - Trln(— Gq 1 + £(*')) + Trh^-G -1 (£')), where ft', 
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Fig. 2. (Color online) Variational parameter /iaf dependence of 
the grand potential for various values of 6 [Eqs. (9) and (14)] for 
U = 1.44 eV, J/U = 0.15, and the realistic set of parameters for 
Sr 2 Ir0 4 listed in Table I. 



E(t'), and G(i') are the ground potential, self-energy, and 
Green's function of the reference system H' , respectively. 
Go is the non-interacting Green's function of H. The 
variational condition dCl(t')/dt' = determines an ap- 
propriate reference system H' which describes the origi- 
nal system H approximately. 

To study the symmetry-broken long-range-ordered 
states in the VCA, we introduce suitably chosen ficti- 
tious Weiss fields in a set of variational parameters t' . 
In this study, we introduce the Weiss field acting on the 
Kramer's doublet 



H AF = h AF J2e lQr 4t 



(14) 



with Q = (tt, 7t). Note that we only consider the Weiss 
field corresponding to the in-plane antiferromagnetic or- 
dering since the energy for the states with the out-of- 
plane antiferromagnetic ordering is always higher than 
that for the in-plane one, as reported previously. 5 ' 

Fig. 2 shows Q,(h,AF,0) — £1(0,0) per site for different 
values of and with U = 1.44 eV, J/U = 0.15, and the 
realistic set of parameters for Sr2lr04 listed in Table I. 
The fact that this quantity has a minimum at a finite 
value of /iaf indicates that the ground state is antiferro- 
magnetically ordered. Carrying out careful calculations, 
we also find that the optimal value of is w 25° for 
this set of model parameters. 

Let us now study the one-particle excitations using 
the Green's function optimized above. The results for 
the momentum resolved one-particle excitation spectra 
and the density of states are shown in Fig. 3. Here, to 
understand the effects of A as well as U , we choose three 
different sets of parameters with (a) and (d): A = U = 
J = 0, (b) and (e): A = 0.37 eV and U = J = 0, and (c) 
and (f): A = 0.37 eV, U = 1.44 eV, and J/U = 0.15. 

As clearly seen in Fig. 3 (a), when the SOC and the 
local Coulomb interactions are absent, the bands con- 
sist of two one-dimension-like narrower bands (ranging 
from — 1 eV to eV) , which are originated from d yz and 
d Z x orbitals, and the remaining band with broader band 
width extending from —3.2 eV to 0.5 eV, which has a 
characteristic of d xy . As shown in Fig. 3 (d), the pro- 
jected density of states onto J cr f = 1/2 extends to the 
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Fig. 3. (Color online) The momentum resolved one-particle excitation spectra (upper panels) and the density of states p{ui) (lower 
panels) for (a) and (d): A = U = J = 0, (b) and (e): A = 0.37 eV and U = J = 0, and (c) and (f): A = 0.37 eV, U = 1.44 eV and 
J/U = 0.15. Other parameters used are the realistic ones for Sr2ir04 listed in Table I. V, M and X correspond to (k x ,k y ) = (0,0), 
(n, 0), and (tt,tt) in first Brillouin zone, respectively. In (d), (e), and (f), black solid (red dashed) lines indicate the total density of 
states (the partial density of states projected onto J e jf = 1/2). u) = corresponds to the Fermi level. Delta functions are represented 
using Lorentzian functions with its half- width at half- maximum r\ = O.lti for (a), (b), and (c), and r\ = 0.005ti for (d), (e), and (f). 



whole energy band, indicating that J e s = 1/2 is not a 
good quantity to describe the one-particle excitations. 
When the SOC is turned on in Fig. 3 (b) and (e), the 
three different bands are observed, i.e., the highest one 
located from —1.1 to 0.6 eV in energy, the intermedi- 
ate one ranging from —1.4 to 0.2 eV, and the lowest one 
extending from —3.2 to —0.2 eV. For this case, we can 
observe that while the partial density of states projected 
onto J ff = 1/2 is still extended to the whole energy re- 
gion, most of the unoccupied states (from ~0.2 to ~0.6 
eV) has a characteristic of J e g = 1/2. Finally, when we 
include the Coulomb interactions in Fig. 3 (c) and (f), we 
can clearly see that the "upper Hubbard band" , which is 
located above the Fermi level, is well separated from the 
valence band. We can also see in Fig. 3 (f) that almost 
all the unoccupied states is of J e g = 1/2 characteris- 
tic, which provides a numerical evidence that the ground 
state is well characterized by the J e g — 1/2 Mott insu- 
lator. This finding is in good qualitative agreement with 
the recent experiments on S^IrOi. 4 ) 

5. Conclusions 

We have employed the exact diagonalization method 
to study the ground state phase diagram for the three- 
band Hubbard model with the SOC. We have found that 
in the Mott insulating phase the ground state can be well 
described by the Kramer's doublet. This suggests that 
no apparent symmetry change exists between the two 
extreme states, i.e., the S =1/2 Mott insulator and the 
J c s = 1/2 Mott insulator, which appear for A = and 
oo, respectively. We have also studied the one-particle 
excitations using the VCA, and found that both A and 



U are essential to realize the J e s = 1/2 Mott insulator, 
where the unoccupied states are mostly of J e g = 1/2 
characteristic. 
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